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We propose new algorithm for numerical modelling of Bose-Einstein correlations (BEC) . It is based 
on the fact that identical particles subjected to BEC do, by definition, bunch themselves in a 
maximal possible way, restricted only by conservation laws, in the same cells in phase-space. 
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The study of Bose-Einstein correlations (BEC) 
has already long history It is motivated mainly 

by their ability to provide space-time information about 
multiparticle production processes. This is of particular 
importance in searches for proper dynamical evolution of 
heavy ion collisions, especially with respect to the pro- 
duction of Quark Gluon Plasma (QGP) Q]. Because of 
their complexity the multiparticle production processes 
are modelled using Monte Carlo event generators of dif- 
ferent sorts Their totally probabilistic structure 
excludes a priori genuine BEC, which are of purely quan- 
tum statistical origin. One can only hope to model BEC 
by changing the original output of these generators in 
such a way as to reproduce measured experimentally two- 
particle correlation function. 
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defined as the ratio of the two-particle distributions to 
the product of the single-particle distributions. This is 
achieved by a suitable bunching finally produced identi- 
cal particles (mostly pions of the same sign) in phase- 
space and can be performed in many ways. In |^] 
it is done locally by shifting slightly (pulling closer to 
each other) momenta of like-charge mesons according to 
some prescribed weights and correcting afterwards for the 
energy-momentum imbalance introduced this way. Other 
method is to select events already showing such bunch- 
ing and count them as many times as necessary in or- 
der to get desirable C2{Q)- This is achieved by using 
special weights calculated for each event from the out- 
come of the event generator used iQ-^ . Such procedure 
preserves energy-momentum balance in each event but 
changes final single-particle distributions obtained from 
the event generator and should be corrected accordingly. 
Because in practical terms it is prohibitively time con- 
suming, sometimes this method is reduced to the simple 
multiplication of particle configuration recorded in each 
event by such weight (to this cathegory belongs also " af- 
terburner" method ||lo| , |ll[| ) . However, in this case energy- 
momentum balance is again altered and should be cor- 
rected accordingly. 

In all these approaches parameters describing appro- 
priate weights do not reflect in any way the possible 



structure of the hadronizing source. They have been 
chosen simply to describe the experimental C2 function 
(|l|). In this letter we would like to propose another 
method of introducing desired bunching of identical par- 
ticles. It will make direct use of the fact that identi- 
cal particles subjected to BEC have, by definition, very 
strong tendency to bunching themselves in a maximal 
possible way (restricted only by conservation laws, es- 
pecially by the energy-momentum conservation) in the 
same cells in phase-space Usually event generator 

provides us in each event with a number of charged and 
neutral secondaries (assumed to be pions, for simplic- 
ity) of which we know their energy-momenta and some- 
times also spacio-temporal positions of their production 
points. The weighting procedure mentioned above is just 
a kind of filtering of only those events which, by shear 
chance, are already showing desired pattern of BEC. Per- 
forming such filtering already in each event would re- 
sult in BEC pattern preserving at the same time both 
the energy-momenta and total multiplicity distributions. 
This is possible if one resignes from some information 
provided by event generator. The energy-momenta and 
total charges cannot be changed because they are directly 
measured. However, the spacio-temporal pattern of the 
event or the charge distribution among particles in the 
event can be altered as not directly observable. Changes 
in spacio-temporal pattern corresponds in a sense to in- 
troduction of quantum mechanical clement of uncertainty 
to the otherwise classical event generator and has been 
discussed to some extend in We shall concen- 

trate therefore in what follows on charge allocation to 
the produced particles. In particular, we shall propose to 
subject all identical particles produced in a given event to 
a kind of specific charge assignement filter, which either 
allocates charges to particles, if event generator does not 
perform this function, or otherwise changes their origi- 
nal charge allocation. This is quite different approach 
from that usually discussed , because instead of sym- 
metrization of appropriate multiparticle wave function 
one works in the number of particles basis and imple- 
ments bosonic character of the produced secondaries by 
bunching the like-ones in the phase space. In it was 
done by introducing concept of elementary cmmitting cell 
groupping particles of the same charge {EEC), in by 
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using a new statistical model based on the information 
theory approach. In any case (following the original HBT 
effect 0]) one has to model quantity measuring correla- 
tions of fluctuations present in the system: 

(rii?i2) = {ni){n2) + {{ni ~ (ni)) (712 - (^2))) 

= {ni){n2) + p(j{ni)a{n2). (2) 

Here cr{n) are dispersions of the multiplicity distribu- 
tion P{n) and p is the correlation coefficient depending 
on the type of the particles produced: p = +1,-1,0 
for bosons, fermions and Boltzmann statistics, respec- 
tively. The proposed algorithm should then provide us 
with C2(Q), which can be written in the form: 



C2{Q \pi -Pj\) 



= 1 
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It means that BEC can be regarded as a measure of cor- 
relation of fluctuations. To get p > it is enough to 
select one of the produced particles, allocate to it some 
charge, and then allocate (in some prescribe way) the 
same charge to as many particles located near it in the 
phase space as possible. In this way one forms a cell 
in phase-space occupied by like-charged particles only. 
This process should be repeated untill all particles are 
used. Notice that, contrary to all previous methods of 
modelling BEC, neither the original energy-momentum 
distributions nor the spacio-temporal pattern provided 
by our event generator are altered. On the other hand, 
we change completely the already existing charge pat- 
tern but retain both the initial charge of the system and 
its total multiplicity distribution. Therefore this method 
works only when we can resign from controling the charge 
flow during hadronization process. The procedure of for- 
mation of such cells will be controlled by appropriate 
weights deciding whether or not a given neighbour of the 
initially selected particle should be counted as its an- 
other member. Using selection procedure which leads to 
a geometrical particle distribution in cells, one maximizes 
second term in (^ because in this case a —< n >. This 
is the general idea of what should be done. 

The proposed algoritm is then as follows: let in the 
p-h gvent our generator provides us with ni = n 



{+) 



particles. Keeping their energy- momenta {pj} 
and spacio-temporal positions {xi \ intact, we allocate to 
them charges in the following way: 

(1) The SIGN is choosen randomly from: 



' or " 0" pool, with weights proportional to p\ 



(+) 



n[^^/ni, p[ ^ 



,(-) 



/ni and p\ 



(0) 



,(0) 



/ni. It is 



attached to particle (i) chosen randomly from the 
particles produced in this event and not yet reas- 
signed new charges. 

(2) Distances in momenta, 6ij (p) — \pi ~ Pj\, between 
the chosen particle («) and all other particles (j) 



still without signs are calculated and arranged in 
ascending order with j = 1 denoting the nearest 
neighbour of particle (z). To each 5ij{p) an appro- 
priate weight P{i,j) G (0, 1) is then assigned (the 
form of which will be discussed below). 



(3) A random number r G (0, 1) is selected from a uni- 
form distribution. If nf^'-^^ > 0, i.e., if there are 
still particles of given SIGN with not reassigned 
charges, one checks the particles (j) in ascending 
order of j and if r < P{i,j) then charge SIGN is 
assigned also to the particle (j), the original mul- 
tiplicity of particles with this SIGN is reduced by 
one, nf^'^^ = nf^'^^ — 1, and the next particle is 
selected: (j) ^ {j + 1). However, if r > P{i,j) or 



^SIGN 

updated values oi p\^\ p[~^ and pf \ 



finishes when 



then one returns to point (1) with the 

Procedure 



= 0, in which 
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case one proceeds to the next event. 

It is easy to check that this algoritm indeed leads to 
a geometrical (Bose-Einstein) distribution of particles 
in the phase-space cells formed in this procedure p^ ] 
accounting therefore for their bosonic character (Bose- 
Einstein statistics). That this leads to C2{Q) > 1 
for Q — > was already demonstrated in statistical 
model based on information theory (with conservation 
of charges imposed) [Q. In contrast to case we are 
working with particles already provided us by some event 
generator and both the sizes of phase-space cells and their 
number are varying from event to event and within given 
event depending on the values of weights P{ij). It is im- 
portant to realize that, because we do not restrict a pri- 
ori the number of particles which can be put in a given 
cell, we are automatically getting BEC of all orders. It 
means that C2{Q = 0) calculated in such environmet of 
multiparticle BEC can exceed 2 (cf., [^). 

To demonstrate abilities of our algoritm we show in 
Fig. 1 results for two different models of hadronization 
of mass M: CAS and MaxEnt (limited for simplicity to 
one-dimensional cases only) . CAS is based on the phase- 
space and space-time cascade model discussed by us re- 
cently |T^, MaxEnt is simple statistical model based on 
information theory concepts providing us only with 
the momenta distributions. None of them shows BEC. As 
can be seen, our algoritm results in a clear BEC type be- 
haviour of correlation functions C2{Q) (which come out 
to be of exponential shape) even with constant weights 
P{i,j) — 0.5. They apparently do not depend much 
on the type of hadronization used. The comparison be- 
tween models is done in the following way: to each CAS 
event characterized by the multiplicity ni one builds the 
corresponding MaxEnt event according to the procedure 
outlined in p^ ] (calculating the corresponding Lagrange 
multipher fii or "temperature" Ti — I /Pi, which describes 
distribution of particles in phase space). Using now the 
same multiplicities, ni, ' and nf\ as in CAS 

one calculates the corresponding BEC in MaxEnt. 
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FIG. 1. Comparison of BEC obtained using our algoritm 
for CAS (left panels) and MaxEnt (right panels) types of 
the emmiting one-dimensional sources of different masses AI. 
Black symbols are for P = 0.5 and open symbols for the "most 
natural" choices of P discussed in text. 

The only additional place where dynamics can enter 
are weights P{ij)- Leaving aside more detailed analysis of 
this problem, we shall display here, in addition to the con- 
stant weigths arbitrary chosen to be equal P{ii) — 0.5, 
also a kind of the "most natural" choices of these weights, 
which uses only the available information provided by 
event generator: 



P{ij) = exp 



(4) 



for CAS (particles (ij) described by spacio-temporal 
wave packets would have widths given by their momen- 
tum separation Sij{p)) and 



P{ij) = exp 



(5) 



for MaxEnt {fix is transverse mass put here to be equal 
0.3 GeV and the role of spatial dimension is now played 
by the "temperature" Ti of the Z*'' event mentioned 
above). The BEC effect is now weaker, but evidently 
present. 

The similarity of BEC patterns in both types of models 
originates in the fact that BEC effect is in our case given 
entirely by the number of particles of the same charge in 



a phase-space cell. This depends on P, the bigger P the 
more particles and bigger C2iQ = 0); smaller P leads to 
the increasing number of cells, which, in turn, results in 
decreasing C2{Q = 0), as already noticed in Only 
after connecting P with details of hadronization process, 
some differences between models start to be visible. They 
originate entirely in differences in the number of elemen- 
tary cells and number of particles located in them. There- 
fore the "sizes" R obtained from the exponential fits to 
results in Fig. 1 (like C2{Q) ~ l-l-A-exp(— Q-i?) where A 
being usually called chaoticity parameter |^) correspond 
to the sizes of respective elementary cells rather than 
to sizes of the hadronizing sources itself. For P = 0.5 
the "size" R varies weakly between 0.66 to 0.87 fm from 
M = 10 to 100 GeV for CAS and between 0.83 and 1.32 
fm for MaxEnt whereas for other weights discussed above 
it varies, respectively, from 0.60 to 0.57 fm for CAS and 
from 0.94 to 1.67 fm MaxEnt. This should be contrasted 
with the real sizes of CAS sources changing from 0.29 fm 
for M = 10 GeV to 1.61 fm for M = 100 GeV. 
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FIG. 2. Dependence of BEC calculated using our algorithm 
for the CAS type of the hadronic source on the actual num- 
ber of sources emmiting particles: the original mass M — 100 
GeV {k = 0) is divided into 2*° sources with k = 1,2,3, re- 
spectively, positioned in the same place. Upper panel is for 
P = 0.5 and lower for the "most natural" choice of P dis- 
cussed in text. 

So far we were considering only single sources. Sup- 
pose now that mass M consists of a number of sources 
hadronizing independently. It turns out that result- 
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ing C2's are sensitive to whether one applies our algo- 
rithm of assigning charges to all particles from subsources [2] 
taken together or to each of the subsource independently. [3] 
Whereas the later case results in the similar "sizes" R 
(defined as before) with C2{Q — 0) = 1 + A faUing 
dramatically with increasing k (roughly like 1/2'' sub- 
sources), the former leads to rougly the same C2{Q = 0) 
but the "size" R is now increasing. The results for this 
case are shown in Fig. 2. The corresponding Rk are now 
equal to, respectivly, 0.87 fm, 1.29 fm, 1.99 fm and 3.35 
fm for P = 0.5 and 0.57 fm, 3.26 fm, 4.01 and 5.59 fm 
for gaussian P as defined above ||l^ . 

To summarize: we propose new and simple method of 
numerical modelling of BEC. It is based on the expected 
phase-space behaviour of identical bosons, conserves the 
energy-momenta and does not alter spacio-temporal pat- 
tern of events or any single particle inclusive distribution. 
We find that value of C2{Q = 0) (defining chaoticity pa- 
rameter A) depends on the number of elementary cells 
in the way already discussed in Q and that "radius" R 
extracted from the exponential fits is practically indepen- 
dent on the size of the source, provided it is a single one. 
In the case when it is composed of a number of elemen- 
tary sources, R increases with their number, unless they 
are treated independently by our algorithm. It is because 
one has in this case higher density of particles what re- 
sults in smaller average Q, and this in turn leads to bigger 
R. On the contrary, for the independently treated sources 
density of particles subjected to our algorithm does not 
change, hence the avarage Q and R remain essentially 
the same. However, because in this case the influence of 
pairs of particles from different subsources increases, the 
effective A = C2(0) — 1 decreases now substantially (as 
was already observed in jlj]). It means that with such 
algoritm we can attempt to fit experimental data. One 
should remember, however, that so far our results are 
limited to one-dimensional case only and most probably 
should be first generalized to 3-dimensional case as well. 
Also the role of resonances (for which charge selection 
should be done independently) and final state interac- 
tions (in particular Coulomb interactions) should be first 
clarified. This will be done elsewhere. Finally, we would 
like to mention that our algorithm leads to strong intcr- 
mittency showing up after its application, in a fashion [15] 
similar to that discussed already in (enhancing the 
natural intermittency pattern existing already in CAS 

and introducing it in MaxEnt case where without [16] 
BEC it was not present). 
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